Relaxation of Pseudo pure states : The Role of Cross-Correlations 
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In Quantum Information Processing by NMR one of the major challenges is relaxation or 
decoherence. Often it is found that the equilibrium mixed state of a spin system is not suitable as 
an initial state for computation and a definite initial state is required to be prepared prior to the 
computation. As these preferred initial states are non-equilibrium states, they are not stationary 
and are destroyed with time as the spin system relaxes toward its equilibrium, introducing error 
in computation. Since it is not possible to cut off the relaxation processes completely, attempts 
are going on to develop alternate strategies like Quantum Error Correction Codes or Noiseless 
Subsystems. Here we study the relaxation behavior of various Pseudo Pure States and analyze the 
role of Cross terms between different relaxation processes, known as Cross-correlation. It is found 
that while cross-correlations accelerate the relaxation of certain pseudo pure states, they retard 
that of others. 



I. INTRODUCTION 

Quantum information processing (QIP) often requires pure state as the initial state . Shor's 
j prime factorizing algorithm jjj, Grover search algorithm ^ are few examples. Creation of pure 
state in NMR is not easy due to small gaps between nuclear magnetic energy levels and demands 
unrealistic experimental conditions like near absolute zero temperature or extremely high magnetic 
field. This problem has been circumvented by creating a Pseudo Pure State (PPS). While in a pure 
state all energy levels except one have zero populations, in a PPS all levels except one have equal 
populations. Since the uniform background populations do not contribute to the NMR signal, such 
a state then mimics a pure state. Several methods of creating PPS have been developed like spatial 
average QQ. logi ca, .abeling BflflQ, te m po r a, averaging Q, spatial averaged logica, la- 
beling technique (SALLT) [12]. However pseudo pure state, as well as pure states are not stationary 
and are destroyed with time as the spin system relaxes toward equilibrium. In QIP there are also 
cases where one or more qubits are initialized to a suitable state at the beginning of the computation 
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and are used as storage or memory qubits at the end of the computation performed on some other 
qubits Q. In these cases it is important for memory qubits to be in the initialized state till the time 
they are in use since deviation from the initial state adds error to the output result. Since it is not 
possible to stop decay of a state which is away from equilibrium, alternate strategies like Quantum 
Error Correction [l^j], Noiseless subspace 15, [r| are being tried. Recently Sarthour et al. [l^ has 



reported a detailed study of relaxation of pseudo pure states and few other states in a quadrupolar 
system. Here we experimentally examine the lifetime of various pseudo pure states in a weakly 
J-coupled two qubit system. We find that cross terms (known as cross-correlation) between different 
pathways of relaxation of a spin can retard the relaxation of certain PPS and accelerate that of others. 

In 1946 Bloch formulated the behavior of populations or longitudinal magnetizations when 
they are perturbed from the equilibrium [18j . The recovery toward equilibrium is exponential for a 



two level system and for a complex system the recover y in vo 
complex systems the von Neumann-Liouville equation {'. 



ves several time constants 19]. For 



2l| describes mathematically the time 



evolution of the density matrix in the magnetic resonance phenomena. For system having more than 
one spin the relaxation is described by a matrix called the relaxation matrix whose elements are 
linear combinations of spectral densities, which in turn are Fourier transforms of time correlation 



function 



22j of the fluctuations of the various interactions responsible for relaxation. There exist 



several different mechanisms for relaxation, such as, time dependent dipole-dipole(DD) interaction, 



221. The 



chemical shift anisotropy(CSA), quadrupolar interaction and spin rotation interaction 
correlation function gives the time correlations between different values of the interactions. The 
final correlation function has two major parts, namely the 'Auto-correlation' part which gives at two 
different times the correlation between the same relaxation interaction and the 'Cross-correlation' 
part which gives the time correlation between two different relaxation interactions. The mathemat- 



ics of cross correlation can be found in detail, in works of Schneider 



Hubbard 



Blicharski 12 5l 



and 



Recently a few models have been suggested to study the decoherence of the quantum 



coherence, the off-diagonal elements in density matrix |27l . |28(. It can be shown that in absence 
of r.f. pulses and under secular approximation the relaxation of the diagonal and the off-diagonal 
elements of the density matrix are independent Here we study the longitudinal relaxation that 
is the relaxation of the diagonal elements of the density matrix and the role of cross-correlations in 
it. 
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II. THEORY 

A. The Pseudo Pure State (PPS) : In terms of magnetization modes 

In terms of magnetization modes the equilibrium density matrix of a two spin system is given by 
B,0,y,Q[FiglH, 

Xeq = llhz + l2hz (1) 

where 71 and 72 are gyro-magnetic ratios of the two spins I\ and h respectively. The density matrix 
of a general state can be written as, 

X = ±a^Jxz ± P^hz ± v~ii2I lz hz} (2) 

which for the condition a / ji=flj2= l/ 'yi=K, corresponds to the density matrix of a PPS given by 

Xpps = K[±I 1Z ± I 2Z ± 2I 1Z I 2Z ] (3) 

where, K is a constant, the value of which depends on the method of creation of PPS. 
The first two terms in the right hand side in Eq|2] and Eq|3] are the single spin order modes for the 
first and second spin respectively while the last term is the two spin order mode of the two spins 
Choosing properly the signs of the modes, the various PPS of a two-qubit system are, 

XZ = K[+I 1Z + hz + 2I lz hz] 
X f ps = K[-I 1Z + hz + 21 lz Iiz\ 
Xp ° ps = K[+I 1Z - hz + 2I lz hz] 

XH = K[+I 1Z + hz ~ 2I lz hz] (4) 

The relative populations of the states for different PPS are shown in Fig. |21 As seen in Eq.2, in PPS 
the coefficients of the all three modes are equal. On the other hand equilibrium density matrix does 
not contain any two spin order mode. To reach EqEl starting from Eq^ the two spin order mode 
has to be created and at the same time the coefficients of all the modes have to be made equal. 

B. Relaxation of Magnetization Modes 



The equation of motion of modes M is given by 



22 



±M(t) = T[M(t) - M(oo)] (5) 
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where f is the relaxation matrix and M(oo) is the equilibrium values of a mode. For a weakly 
coupled two-spin system relaxing via mutual dipolar interaction and the CSA relaxation, the two 
dominant mechanism of relaxation of spin half nuclei in liquid state, the above equation takes the 
form, 



hzit) 




Pi &12 &L.12 




hz(0) - Jiz(oo) 




hzit) 




Cl2 P2 #2,12 




hz{0) - hz{oo) 


(6) 


Ihzhzit) 




01,12 02,12 Pl2 




2hzhz 





where is the self relaxation rate of the single spin order mode of spin i, p^ is the self relaxation 
rate of the two spin order mode of spin i and j, Cy is the cross- relaxation (Nuclear Overhouser Effect, 
NOE) rate between spins i and j and c^y is the cross-correlation term between CSA relaxation of 
spin i and the dipolar relaxation between the spins i and j. p and a involve only the auto-correlation 
terms and 5 involves only the cross-correlation terms |22|. Magnetization modes of one order relaxes 
to other orders through cross-correlation and in absence of it the relaxation matrix becomes block 
diagonal within each order. The relaxation of modes are in general dominated by their self relaxation 
p, but in case of samples having long T\, the cross-correlation terms become comparable with self- 
relaxation and play an important role in relaxation of the spins. 
The formal solution of Eq. El is given by, 

M(t) = M(oo) + [M (0) - M(oo)}exp(-ft) (7) 

As time evolution of various modes are coupled, a general solution of the above equation requires 
diagonalization of the relaxation matrix. However, in the initial rate approximation Eq.7 can be 
written (for small values of t=r) as, 

M(r) = M(oo) + [M(0)-M(oo)][l-fr] 

= M(0) - f r[M(0) - Af(oo)] (8) 

This equation asserts that in the initial rate approximation (for low r), the decay or growth of a 
mode is linear with time and the initial slope is proportional to the corresponding relaxation matrix 
element. If the modes are allowed to relax for a longer time, their decay or growth deviates from 
the linear nature and adopts a multi-exponential behavior to finally reach the equilibrium [22]. 
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C. Relaxation of Pseudo pure state 

Let a two qubit system be in 1 00) PPS at t=0. 

X°°(0) = K[I 1Z + I 2Z + 2I 1Z I 2Z ] (9) 

After time t it will relax to, 

x 00 (t) = (K + Ai(t))iiz + (K + A 2 (t))I 2Z + (K + A 12 (t))2I lz I 2Z (10) 

where Ai(t),A 2 (t) and A 12 (£) are the time dependent deviations of respective modes from their 
initial values. The deviation of the two spin order can be measured from spectrum of either spin. 
EqJTUl can also be written as, 

= (K + A 12 {t))[I lz + I 2Z + 2I 1Z I 2Z ] + (A 1 (t) - A 12 (t))I lz + (A 2 (t) - A l2 {t))I 2Z (11) 

The first term is the pseudo pure state with the coefficient decreasing in time while the other two 
terms are the excesses of the single spin order modes with coefficients increasing in time. For other 
pseudo pure states EqlTTl becomes. 

x °i(t) = (K + A 12 (t))[-I lz + I 2Z + 2I 1Z I 2Z ] + (A 1 (t) + A 12 (t))/ lz + (A 2 (t) - A 12 (t))I 2Z (12) 
x w (t) = (A: + A 12 (t))[/ lz -/ 2Z + 2/iz/2z] + (A 1 (t)-A 12 (t))/ lz + (A 2 (t)+Ai 2 (t))/ 2Z (13) 
X n (t) — (K - A 12 (t))[I lz + I 2Z - 2I 1Z I 2Z ] + (Ai(t) + A l2 {t))I lz + (A 2 (t) + A l2 {t))I 2Z (14) 



In the initial rate approximation (using EqJHJ) we obtain for the 1 00) pps, 

A x (r) = r[ Pl ( 7l - K) + a 12 ( l2 - K) - K8 x , l2 ] (15) 

A 2 (r) = r[ ( r 12 ( 7l - K) + Pl (j 2 - K) - K5 2 , 12 ] (16) 

A 12 (r) = r[5i,i 2 (7i - K) + 5 2 , 12 ( 72 - K) - Kp 12 ] (17) 



Let the coefficients of the PPS term and the two single spin order modes l\ Z and I 2Z in Eq^2 be 
called as A,B and C respectively. FigOU schematically shows the time evolution of the coefficients 
A,B and C for 1 00) PPS. Any coefficient for any PPS at any instant, is simply the initial value plus 
the total deviation due to the auto and the cross-correlations. For example, A for 1 00) PPS at time 
r, is A 00 (t) = K + A° a l to (r) + A™(t), where K is the initial value, and A™ uto (r) and A™(t) are the 
deviations at r due to auto-correlation and cross-correlation parts respectively. 
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(a) Contribution of auto- correlation terms to the deviation 

Putting the values of the deviations of different modes obtained from Eq. (jl5HlTj) in Eq. ([11111%}) . 
we obtain the contribution only of auto-correlation terms to the deviation from initial value of the 
coefficients A,B and C under initial rate approximation (at t=r) as, 

AZtoir) = AZ to (r) = AZ (t) = A^r) = -Kp 12 r 
Kltoir) = [Piili ~K) + a 12 ( 72 - K) + Kp l2 ]r ; B° a l to (r) = [ Pl ( 7 i + K ) + n 12 { l2 - K) - Kp 12 ]r 
BZo(t) = M7i ~K) + a 12 ( 72 + K) + K Pi2 ]t ; B^ io (r) = [ Pl ( 7l - K) + a l2 { l2 - K) + Kp 12 ]r 
CZo(r) = [0-12(71 - K) + p 2 ( 72 - K) + K Pl2 ]r ; Cf uto [ T ) = [<n 2 ( 7 i + K) + p 2 ( 72 - K) + Kp 12 ]r 
CZo(r) = [oia(7i - K) + p 2 ( 7 2 +K) - Kp 12 ]r ; C^ to (r) = [<r 12 ( 7l - X) + p 2 ( 72 - X) + Kp 12 ]r (18) 

It is evident that in absence of cross-correlations the 1 00) and |11) PPS relax at the same initial rate 
since A™ ut o = A a i to , B° a ° uto = B x ^ uto and C°° io = CH to . However the same is not true for |01) and |10) 
PPS. 



(b) Contribution of cross-correlation terms to the deviation 



The contribution by the cross-correlation terms is given by, 

•4LV) = +[*L,ia(7i ~K) + #2,12(72 - K)\r 
A° c l(r) = +Pi,i2(7i + K) + #2,12(72 ~ K))t 
Al° c (r) = +[5i,i 2 (7i -K) + #2,12(72 + K))t 
AIKt) = -0i,i 2 (7i -K) + #2,12(72 - K)]r 

B°cc(t) = -0i,i27i + #2,12(72 - K)\t ; fi£(r) = +0i,i 2 7i + #2,12(72 - K)\r 

Blc(r) = -0i,i27i + #2,12(72 + K)\t ; B£(t) = +0i,i 2 7i + #2,12(72 - 

CV) = -0i,i2(7i -K) + 8 2>12 j 2 ]t ; C£(r) = -0i,i 2 (7i + K) + #2,1272]^ 

C(^) = +0i,i2(7i - + #2,1272]^ ; C^ 1 (r) = +0 ljl2 (7i - K) + 5 2A2l2 ]r (19) 

The important thing is that the presence of cross-correlation can lead to differential relaxation 
of all PPS. Positive cross-correlation rates #1,12 and #2,125 slow down the relaxation of all the 
three coefficients for |00) PPS since (A 00 (r) > ^(r), B°°(r) < BZo(r)X°°(r) < C to (r)), 
while make the relaxation of all three coefficients faster for |11) PPS since (A u (t) < 
Ali t0 (r),B n (r) > JB^Jt), C u (t) > C£Jt)). For |01) and |10) PPS cross-correlations 



give a mixed effect since (A°\r) > A^ {t),B»\t) > ^ to (r), C 01 (r) < C to (r)) and 
(A 10 (r) > AZo(t),B 10 (t) < BZ o (r),C 10 (r) > C to (r)). As the contributions of the anto- 
correlation part for 1 00) and |11) PPS are equal, we have monitored the relaxation behavior only of 
1 00) and |11) PPS to study the effect of cross-correlations. 



For samples having long T 1; where the cross-correlations becomes comparable with auto-correlation 
rates, the four PPS relax with four different rates and the difference increases with the increased 
value of the cross-correlation terms. The three coefficients A,B and C (normalized to the equilibrium 
line intensities) in terms of proton and fluorine line intensities for 1 00) PPS are, 

H°(t) - H\t) F°(t) - F\t) 



A{t) 



W-(pa) +H°(oo) F 1 (oo) + F°(oo) 
2H l (t) 



B(t) 



H°(oo) + H l (oo) 



and for |11) PPS are, 



H\t)-H°(t) F\t)-F°(t) 



ff x (oo) +H°(oo) F 1 (oo) + F°(oo) 
2H°(t) 



Bit) 



H°(oo) + ^(00) 



2F°(t) 

[t) ~ F°(oo) + Fi(oo) (21) 
where, H° and H l are intensities of the two proton transitions, when the fluorine spin is respectively 
in state |0) and |1). Similarly F° and F 1 are intensities of two fluorine transitions corresponding to 
the proton spin being respectively in the state |0) and |1), as shown in Fig^and FigEl H°(t) and 
H°{oo) give the H° line intensity respectively at time t and at equilibrium. Thus by monitoring the 
intensities of the two proton and two fluorine transitions as a function of time, one can calculate the 
coefficient A(t) which is a measure of decay of PPS. 



s 



III. SIMULATION 



Relaxation of the coefficients A,B and C have been simulated using MATLAB, for a weakly coupled 
19 F- l H system. The relaxation matrix used for the simulation is, 

0.3125 0.02 <5 1)12 
0.02 0.33 82,12 

#1,12 #2,12 0.33 

FigH] shows the decay of coefficient A with time, A 00 and .4. 11 show no difference in decay rate in 
absence of cross-correlation rates. As 5^12 and #2,12 are increased more and more difference in decay 
rate is observed. FigJS] shows growth of coefficients B and C. As #2,12 is taken smaller than #1,12, 
difference in decay rate between C 00 and C 11 is found to be less than between B 00 and B n . 



IV. EXPERIMENTAL 



All the relaxation measurement were performed on a two qubit sample formed by one fluorine 
and one proton of 5-fluro 1,3-dimethyl uracil yielding an AX spin system with a J-coupling of 5.8 
Hz. Longitudinal relaxation time constants for 19 F and 1 H are 6 and 7.2 Sec respectively at room 
temperature (300K). All the experiments were performed in a Bruker DRX 500 MHz spectrometer 
where the resonance frequencies for 19 F and 1 H are 470.59 MHz and 500.13 MHz respectively. The 
Pseudo-pure state was prepared by spatial averaging method using J-evolution f|. Relaxation of all 
the three coefficients for 1 00) and |11) PPS has been calculated. Since auto-correlations contribute 
equally to the relaxation of these two PPS, any difference in relaxation rate can be attributed to 
cross-correlation rates. 



Sample temperature was varied to change the correlation time and hence the cross-correlation 
rate S. Four different sample temperatures, 300K, 283K, 263K and 253K were used. Fig JU] shows the 
proton and fluorine spectra obtained using recovery measurement at four different temperatures. 
The spectra correspond to the initial PPS state and that after an interval of 2.5 sec. FigJ7|shows the 
longitudinal relaxation times (Ti) of fluorine and proton as function of temperature obtained from 
initial part of inversion-recovery experiment. A steady decrease in T\ with decreasing temperature 
indicates that the dynamics of the sample molecule is in the short correlation time limit [21]. In 
this limit auto as well as cross- correlations increase linearly with decreasing temperature. 
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All the spectra were fitted to bi-Lorentzian lines in MATLAB and various parameters were ex- 
tracted using the Origin software. Fig|H]shows the decay of the coefficient A calculated independently 
from proton and fluorine spectra. At 300K, A 00 and A 11 showed almost same rate of decay. As the 
temperature was gradually lowered, a steady increase in difference in decay rate was observed. This 
is due to the steady increase in cross-correlation rates with decreasing temperature, which is expected 
in the short correlation time limit. In FigElthe growths of the coefficients B and C are shown. Similar 
to the coefficient A, coefficients B and C also show differences in decay rate between 1 00) and |11) 
PPS at lower temperatures. The difference between B 00 and B n at any temperature was found to 
be larger compared to between C 00 and C n . This is expected since, according to EqE3the dominant 
cross-correlation factor in B 00 and B 11 is 8^x2 which is the cross-correlation between CSA of fluorine 
with fluorine-proton dipolar interaction whereas in C 00 and C n the dominant factor is #2,12 which 
is cross-correlation between CSA of proton, which is much less than fluorine, with fluorine-proton 
dipolar interaction. Thus it is found that at lower temperatures the 1 00) PPS decays slower than 
the 1 11) PPS. The dominant difference in the decay rates arises from the cross-correlations between 
the CSA of the fluorine and the dipolar interaction between the fluorine and the proton spin. To 
the best of our knowledge this is the first study of its kind where the differential decay of the PPS 
has been attributed to cross- correlations. 



V. CONCLUSION 



We have demonstrated here that in samples having long 7\ cross-correlations plays an important 
role in determining the rate of relaxation of pseudo pure state. In QIP sometimes one or more qubits 
having comparatively longer longitudinal relaxation are used as storage or memory qubits. Recently 
Levitt et al. have demonstrated a long living antisymmetric state arrived by shifting the sample 
from high to very low magnetic field, suggesting that this long living state could be used as memory 



qubit |3ll . . In such cases fidelity of computation depends on how much the memory qubits have 
been deviated from the initialized state at the beginning of the computation till the time they are 
actually used. Theoretically it is shown here that in presence of cross-correlations, all the four PPS 
relax with different initial rates. For positive cross-correlations the 1 00) PPS relaxes significantly 
slower than |11) PPS. It is therefore important to choose a proper initial pseudo pure state according 
to the sample. 
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FIGURE CAPTIONS 
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Figure 1. (a) Chemical structure of 5-fluro 1,3-dimethyl uracil. The fluorine and the proton spins 
(shown by circles) are used as the two qubits I\ and I2 respectively, (b) The energy level diagram 
of a two qubit system identifying the four states 00,01,10 and 11. Under high temperature and high 



field approximation 



2l| the relative equilibrium deviation populations are indicated in the bracket 



for each level. Assuming this to be a weakly coupled two spin system the deviation populations 
become proportional to the gyromagnetic ratios 71 and 72. Ij refers to the transition of the j th spin 
when the other spin is in state \k). Thus H° means the proton transition when the fluorine is in 
state |0). 



Figure 2. Population distribution of different energy levels of a two spin system in different 
pseudo-pure states. K is a constant whose value depends on the protocol used for the preparation 
of PPS. (a),(b),(c) and (d) show respectively the |00),|01),|10) and |11) PPS. 

Figure 3. Schematic representation of decay of the coefficient A and growth of the coeffi- 
cients B and C. The magnetization modes are normalized to their respective equilibrium values. In 
each sub-figure the three bars correspond to the modes hz,hz and 2Iizhz from left to right. The 
amount of any mode present at any time is directly proportional to the height of the corresponding 
bar. The numbers provided in the rightmost column represent typical values of the modes, (a) 
Thermal equilibrium. At thermal equilibrium only l\z and I2Z exist, (b) 1 00) pseudo-pure state 
just after creation, where all the three modes are equal in magnitude. For 1 00) PPS all modes are of 
same sign but this is not the case for other PPS [Eq.3]. Coefficient A is the common equal amount 
of all the modes and it is maximum at t=0. (c) The amount of magnetization modes (schematic) at 
time r, after preparation of the PPS at t=0. The two single spin order modes increase and the two 
spin order mode decreases from their initial values, (d) The state of various modes at time r, (same 
as fig.c) redrawn with filled bar to indicate the residual value of A. All the three coefficients A,B 
and C are shown. A (shown by the filled bar), which is the measure of the PPS, has come down by 
the same amount as the two spin order. B (shown by the empty bar) and C (shown by the striped 
bar) are the residual part of the single spin order modes l\z and I2Z respectively, (e) The values of 
various modes and coefficients after a delay r' > r. 



Figure 4. Simulation of decay of coefficient A. The boxes (□) and circles (o) correspond to 
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the |11) and 1 00) PPS respectively. In each plot deviation from initial value (A(t) — „4(0))has been 
plotted. 

Figure 5. Simulation of growth of coefficient B and C. The boxes (□) and circles (o) corre- 
spond to the 1 11) and 1 00) PPS respectively. 

Figure 6. Relaxation of Pseudo pure state as monitored on (a) fluorine spin and (b) proton spin of 
the 5-fluro 1,3-dimethyl uracil at four different temperatures. The top row in (a) and (b) show the 
equilibrium spectrum at each temperature. With decrease in temperatures the lines broaden due 
to decreased T 2 . The second row in (a) and (b) show the spectra corresponding to the 1 00) PPS, 
prepared by spatial averaging method using J-evolution. The state of PPS was measured by 90° 
pulse at each spin. The third row in (a) and (b) show the spectra after an interval of 2.5 seconds 
after creation of the 1 00) PPS. The fourth row shows the spectra immediately after creation of |11) 
PPS and the fifth row, the spectra after 2.5 seconds. 

Figure 7. Longitudinal relaxation time 7\ of fluorine (a) and proton (b) as function of tem- 
perature, measured from the initial part of inversion recovery experiment for each spin. 

Figure 8. The deviation from initial value (at t=0) of the coefficient A of the PPS term cal- 
culated from Proton (left column) and Fluorine (right column) at four different sample temperature. 
The empty (o) and filled (•) circles correspond to the 1 00) and |11) PPS respectively. 

Figure 9. The growth of the coefficients B and C at different sample temperatures. B was 
calculated from Fluorine spectrum while C was calculated from the Proton spectrum. The empty 
(o) and filled (•) circles correspond to the 1 00) and |11) PPS respectively. 
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